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We study the quasiclassical magnetotransport of non-interacting fermions in two dimensions mov- 
ing in a random array of strong scatterers (antidots, impurities or defects) on the background of a 
smooth random potential. We demonstrate that the combination of the two types of disorder in- 
duces a novel mechanism leading to a strong negative magnetoresistance, followed by the saturation 
of the magnetoresistivity pxx{B) at a value determined solely by the smooth disorder. Experimental 
relevance to the transport in semiconductor heterostructures is discussed. 
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Magnetotransport in a two-dimensional electron gas 
(2DEG) has been the subject of intensive research dur- 
ing the last two decades. This interest has been moti- 
vated by the progress in preparing high-quality semicon- 
ductor heterostructures, opening up new areas in both 
fundamental physics and applications, see for a re- 
view. Within the quasiclassical approach (valid for not 
too strong magnetic fields B), impurity scattering is com- 
monly described by a collision integral in the Boltzmann 
equation. This leads, for an isotropic system, to the B- 
independent Drude value of the longitudinal resistivity, 
Pxx{B) = po = m/e'^neT, where rie is the carrier density, 
m the effective mass, and r the transport scattering time. 

It has become clear, however, that this description is 
not always valid. In particular, in the case of smooth 
disorder memory effects induce a strong positive magne- 
toresistance (MR) 1^ followed by an exponential fall-off 
of PxxiB) due to adiabatic localization of drifting elec- 
trons To our knowledge, these effects have not 
been experimentally observed in the electron transport 
in low magnetic fields, since the Shubnikov-dc Haas os- 
cillations develop at lower magnetic fields. On the other 
hand, the memory effects do show up in the composite 
fermion transport, explaining the peculiar shape of the 
MR around half-filling of the lowest Landau level |^,||. 

In the present paper we study the quasiclassical MR of 
a 2DEG moving in a random array of rare strong scat- 
terers (modeled by hard disks) and subject additionally 
to a smooth random potential. Apart from the purely 
theoretical interest, our work has been motivated by two 
types of experimental realizations of this problem. The 
first one is random antidot arrays. Experiments on this 
kind of structures show a strong negative MR 

which has not been analyzed theoretically. Less obvi- 
ously, our model is relevant to transport in the unstruc- 
tured high-mobility 2DEG. To clarify this point, we recall 
that in order to increase the 2DEG mobility the donors in 
currently fabricated heterostructures are separated by a 
large distance kp^ (with kp the Fermi wave number) 



from the 2DEG plane. It is usually assumed that these 
remote donors constitute the main source of disorder, in- 
ducing small-angle scattering of electrons. It is known, 
however, that in samples with a wide spacer (d > 70 nm) 
large-angle scattering on residual impurities |pd|-p^ and 
interface roughness becomes important, limiting the 
mobility with further increasing width of the spacer. In 
particular, Ref. jl^ concludes that about 50% of the re- 
sistivity is determined by such scattering processes, while 
Ref. finds that for samples with very high mobility 
this value is as high as 90%. 

We thus consider the following two-component model 
of disorder: (i) randomly distributed hard-core scatter- 
ers (which we will term "antidots" or "impurities" below) 
with density ns and radius a (where ng^^"^ ^ a 3> kp^), 
and (ii) smooth random potential (correlation radius d, 
momentum relaxation rate tJ^^ , transport mean free path 

= vptl). The mean free path for the scattering on an- 
tidots is l^g^ = vpT^s^ = l/2nsa, while the correspond- 
ing transport mean free path may be somewhat different. 
Is = vpTs = 7^5* with 7 ~ 1 (7 = 3/4 in the model of 
specularly reflecting disks). We will set 7 = 1 for quali- 
tative estimates. We will further assume that tl ^ T5, 
so that the zeio-B resistivity po is determined by the 
hard scatterers, r^^ = r^^ +'''s^ — '^s^ ■ Finally, we will 
assume that the motion in the smooth disorder is not 
adiabatic, i.e. has the form of the guiding center diffu- 
sion (rather than drift). The condition for this is (5 ^ d, 
where 6 is the guiding center shift after one cyclotron 
revolution [see Eq. (||)] . For currently fabricated samples 
this assumption is usually satisfied in the whole range of 
applicability of the quasiclassical theory. 

We start the analysis of the problem by recalling the 
results for the case when only hard scatterers are 
present (r^, = 00), known as the Lorentz gas. In this 
limit the resistivity Pxx{B) reads 

Pxx{B)Ip^ = (1 - e-^^'-^^s ')T{uj,rs), (1) 
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where lUc = eB/mc is the cyclotron frequency and 
J-{x) is a function of order unity with the asymptotics 
< 1) = 1 and T{x > 1) = 7. The first factor 
in Eq. (|^) is nothing else but the fraction of particles 
moving in rosette-like trajectories around the impurities 
and hitting a new impurity with the mean free time t^''^ . 
The rest of the particles do not hit scatterers at all. In 
the sequel, we will only consider classically strong mag- 
netic fields, ujc ^ lOt = 2'K/Tg'\ where the resistivity 

(1) shows a l/B fall-off. Equation ^ is valid below the 

1/2 

percolation threshold, ujf. < uj-p^ic — i-QTvprig (note 
that tOpcrcTs ~ ng^^'^a~^ ^ 1). For larger magnetic 
fields, LUc > Wpcrc, the resistivity is exactly zero, since the 
rosette-like families of cyclotron orbits associated with 
each impurity fail to form an infinite cluster. 

Clearly, adding the long-range disorder will increase 
the diffusion constant D^^ and thus, in the limit lUcTs 3> 
1, the longitudinal resistivity, by setting free those par- 
ticles which are localized in cyclotron orbits not hitting 
impurities. We will be interested in the case of a suffi- 
ciently strong smooth disorder modifying the result (|^) 
in an essential way. Specifically, we will see below that 
new physics emerges in the regime (5 3> a, where 



(2) 



Let us first outline this new physics on a qualitative 
level. Naively, one could think that for S ^ a the resistiv- 
ity should take its Drude value. Indeed, let us associate 
with the particle trajectory a strip of width 2a surround- 
ing it. The particle will hit an impurity if the center of 
the latter is located within this strip. Clearly, in one cy- 
clotron revolution the particle "explores" in this way the 
area 2vFa x 2'k/uJc, the same as it would explore in the 
same time for i? — 0. For 6 ^ a the area explored in 
the second revolution will overlap only weakly with that 
explored in the first one, so that one could think that 
the exploration rate is essentially the same as for _B = 0, 
leading to the mean time ~ r^°^ between the collisions, 
and thus to pxx — Pq- This consideration is, however, 
incorrect, since it neglects memory effects. Specifically, 
there is a probability Pi ^ a/S that the strip after the 
first revolution covers again the starting point (or, in 
other words, there is typically a small relative overlap 
^ a/S). In view of the diffusive dynamics of the guid- 
ing center, the return probability decreases rather slowly 
with the number of revolutions, P„ — Pi/^/n. There- 
fore, the total return probability P — ^^^i Pn is deter- 
mined by the upper cut-off iV, so that the memory effect 
is much stronger than one might at first sight expect, 
(a/(5)iVi/2 > a/5. 

When the memory effect leads only to a small cor- 
rection to the Drude value (the corresponding condition 
will be specified below), the upper cut-off is given by the 
number of cyclotron revolutions it takes for the particle 
to hit the next impurity, N = ujcTg''^ /2tt, so that the 



total return probability is P '--^ {a/S){u!cTsy^'^ ■ This de- 
termines the fraction of the area explored twice, implying 
an effective reduction of the exploration rate and thus a 
negative correction to the resistivity. 



Apxx/po ^ -{a/5){ujcTsY^'^ ^ -(wc/wo)^. 



(3) 



where tJo ~ vpia^lsh) ^^'^ ^ ^^pciSs /hY^^ < ^^porc- 

We turn now to a more rigorous and quantitative 
derivation. Generalizing the formalism of |^ to the 
present case, we start from the following Liouville- 
Boltzmann equation for the distribution function g{Y, <j)) 
of electrons on the Fermi surface, 



(Lo + SL)g{r, (p) cos{(t> - (j)E), 

d I d'^ 

Lq = vpnV + ujc^ ^-2, 

0(p TL ocp'^ 



(4) 
(5) 
(6) 



where n — (cos 0, sin 0) is the unit vector determining 
the direction of velocity, v = vpn, and 0b is the po- 
lar angle of the electric field. The operator Lq describes 
the motion in the smooth random potential, while SL 
corresponds to the scattering on antidots with (random) 
positions R^. The explicit form of the collision operator 
/r (for a hard- wall scatterer) can be found in , but we 
will not need it. We will only use the following properties 
of /r: (i) it is non-zero only within the distance a from 
the impurity location R, and (ii) its Fourier-transform 
/q satisfies /on — (— 1/715x5)11. 

Expanding in SL, averaging over the positions Ri 
of scatterers, and resumming the series, one finds 
the Green's function of the Liouville operator, {{Lq + 
5L)-^) = {Lq + M)-\ where M is the "self-energy" op- 
erator. The resistivity is expressed in terms of M as 
follows (I, 



{m/e^ne){T^^ +Mxx), 



(7) 



where Mxx — J {dip/n) cos (j> M cos (j>. One may use a di- 
agrammatic technique analogous to that developed for 
the Lorentz gas |l4|JT5|] to calculate M. The leading 
term corresponds to a single scattering process, yield- 
ing M^^^ = —nslo and thus the Drude contribution, 
Mx]^ = Tg^ , to the resistivity (Q). The next-order contri- 
bution representing the correction induced by the return 
process is obtained as 



XX 



-ns I — cos( 

TT 



IrD^^^ Ircosc 



(8) 



where D^-^^ — {Lq + Af'^^)^^ is the electron propagator 
with the leading-order self-energy included. 

Since (5 3> a, the propagator D'^^^ (r — r', 0, 0') describ- 
ing propagation from the point r',(/)' to the point r, 
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can be replaced in (||) by D^^\0, cj), cj)'). Furthermore, we 
note that once the particle hits an impurity, its guiding 
center is shifted by an amount of order of the cyclotron 
radius Re- As a result, the contribution of such trajec- 
tories to 

D(i) can be neglected and only non-colliding 
orbits should be taken into account. Since the motion 
without collisions is limited by the times ~ T5 <C tl, 
the particle will return with almost the same direction 
of velocity, i.e. I?^^^ (0, 0, 0') will be peaked at (/) « (j)' . 
We can thus approximate D^^\0, 0, 0') by £)q^^(5(0— 0'), 
where d'^^P — J d(j)D^^^ (0, cf), ((>'). This quantity is easily 
found to be 

p(i) f. e-W^-4\ (a;,.("))V^ 

n—l 

Substituting this in (p|) and using that /on = 
-{l/nsTs)n, wc find = -(I/tistD^^^^ This im- 

plies, according to (0) , the negative MR, 

^Pxx/Po = MH^/Ts^ = -CjI/ujI IVc < Wo, (10) 

where 

u;o = i2nns)'^^VFi2-fls/lLy/\ (11) 

in agreement with the above qualitative considerations. 

This derivation is valid as long as the correction re- 
mains small, i.e. for lOc <C i^o- For stronger magnetic 
fields, when the quantity P defined above becomes large, 
P 3> 1, it acquires the meaning of the number of returns. 
It should then be found self-consistently. Specifically, the 
sum over n is now cut off at the time r'g ~ Pts, since mul- 
tiple returns to the same area lead to a corresponding in- 
crease of the time needed to hit a new impurity. We thus 
get a self-consistency equation, P ^ {a/5){ujcPTsy^^ , 
yielding a 1/ drop of the resistivity, 

Pxx/PO ^ Ts/t's 1/P {Ujo/Ujcy, Wc>WO- (l2) 

It remains to analyze the conditions of validity of 
Eqs. (^^, (p^. First of all, we assumed that the par- 
ticle finds a new scatterer by exploring the new area in 
the course of the diffusive motion of the guiding cen- 
ter. There exists, however, a competing mechanism, 
namely that of the rosette- like motion around a scatterer, 
which determines the transport in the pure Lorentz gas 
{II = 00). Comparing ( [l^ ) with the Lorentz-gas result 
Pxx/pq = jiOcTs [see Eq. (Q)], we find that the two 
formulas match at (5 ~ a. The corresponding crossover 
frequency is i^aoss — ''^F(47rn|.Z|Z£^)"'^/'^. Secondly, we as- 
sumed nsRc ^ 1, or, equivalently, uJc ^ ijJpcrc- It is easy 
to see that in the opposite limit, nsR^ ^ 1, the resis- 
tivity will be determined by the smooth disorder, with 
scattering on antidots giving a small [~ (wporc/wc)^] cor- 
rection only, so that Pxx{B) will have a plateau with the 

value Pxxi^c > Wperc) = m/e^rieTL. 




FIG. 1. Magnetoresistivity at fixed rs and different tl- 
curve A - Lorentz gas (tl = 00), curves B and C with 
rj^^^ > T^"' correspond to regimes B and C, respectively (see 
text). The dotted lines denote the asymptotics (p^. 

Comparing the characteristic frequencies, Ur, wo, 
Across, and Wpcrc, WC conclude that the following three sit- 
uations can be distinguished, depending on the strength 
of the smooth disorder (Fig. |l|): 

(A) uq <C uJr, or, equivalently, Il/Is ^ 
(l/27r^)(ri5^|)^. In this case the smooth disorder hardly 
affects the Lorentz-gas result (|l|). 

(B) < cjo < t^cross < Wpcrc, Or, cquivalcntly, 
2.7(715^1)1/2 ^ i^/ig ^ (^/27r2)(ns;|)2. This is an in- 
termediate situation; the resistivity drops first according 
to Eqs. (|lO|), ( p^ ) and then crosses over at coc Wcross to 
the Lorentz-gas behavior (|^). 

(C) cj^ < Wo < Wporc < Across, Or, cquivalcntly, 
h/h < 2.7{nsliy/^. In this case the Lorentz-gas be- 
havior (P is completely destroyed, and the results (p^), 
( p^ ) hold in the whole range of Uc below Wpcrc- 

For coc > Wporc the resistivity shows in all the cases 
a plateau, as explained above. On the side of strong 
magnetic fields this plateau will be modified either by 
entering into the adiabatic regime at i5 '--^ d, meaning 
Wc ~ Wad = vf{^t^/<PIlY^^ (its implications for the 
present problem will be considered elsewhere jl^) or by 
the development of Shubnikov-de Haas oscillations. 

It is worth mentioning that while we used the con- 
dition Tl 3> T5 for the derivation of our main results, 
a pronounced negative MR will also be observed for 
Tl ^ Ts (which seems to be frequently the relevant situa- 
tion for high-mobility structures (ij]). In this case there 
is a crossover from pxx{0) — ("^/e^'^e)('''L ^ + ^s^) 
Pxxi^c 3> Wpcrc) — in/e^rieTL which takes place around 
Wc ^ Wpcrc (note that wq ^ Wpcrc for tl ~ ts). 

We have performed numerical simulations of the MR 
by fixing parameters of the Lorentz gas (wpcrc/wr — 5.3) 
and the correlation length of the smooth disorder {d/a ~ 
2.5) and changing the strength of the latter. The results 
shown in Fig. ^ are in good agreement with the analytical 
predictions. It is seen that a very weak smooth disorder 
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(giving negligible contribution to po) affects crucially the 
MR. Some deviations in the values of loq from Eq. ( |ll| ) 
(see inset) can be attributed to the fact that the condi- 
tions lUt i^o <^ Wporc^ad are fulfilled only marginally 
in our simulations. 




FIG. 2. Magnetoresistivity at fixed rg and difi^erent tl from 
numerical simulations; tl/ts = oo (Lorentz gas, A), 111 (O), 
70 (□), 37 (O). Inset: wo determined from the fit to Eq. @; 
the full line corresponds to the analytical result (^ij). 

Finally, let us estimate the characteristic values of 
B for existing experiments. Typical parameters in the 
experiments Q on the antidot arrays were rig = 5 x 
10^^ cm~^, ns — (0.6 ^m)^^. Is = 1.3 /im, II = 16 /im. 
This implies the following values for the characteristic 
magnetic fields (we use the obvious notations Br = 
{mc/e)LUr etc.): Br ~ 0.5 T, Bq ~ Bperc ^ 0.3 T. We 
see that the condition of a dilute antidot array assumed 
above, Br <C Bpcrc, is not met. Clearly, in this situation 
the above analytical formulas are not valid quantitatively. 
On the qualitative level, we can conclude that there 
should be a strong fall-off of pxx around B sa 0.3 0.5 T, 
in agreement with experimental results A similar 
negative MR was observed in other experiments with dis- 
ordered antidot arrays l^J^,[l0|. 

We turn now to unstructured high-mobility samples. 
Using the parameters of |l^ {rie — 2 x 10^^ cm^^, mo- 
bility p = 10^ cm^/V • s), we find the mean free path 
I ~ 80pm. Let us assume, following the conclusion of | |T^ ] 
that the zero-field mobility is determined by background 
impurities {i.e. I ~ Is), while II — Wis- Using the typi- 
cal volume concentration of the residual impurities p2[ , 
n^^'' ~ 2.5 X 10^ cm^'^, and the value of the Bohr radius 
in GaAs, ~ 10 nm, we estimate the sheet density of 
strong scatterers as ns ~ asn^s^^ — P'm)~^- With 
these parameters, we find Br — 5 mT, Bq ~ -Bporc — 
60 mT. The condition of the diluted array of scatterers, 

Br<^B 

pore J is now well satisfied. Since Bq ~ Bperci the 
Lorentz-gas behavior (|^) is fully destroyed (the case (C) 
in our classification), and the negative MR is determined 
by the interplay of smooth disorder and strong scatter- 



ers, as described above. Note that though parametri- 
cally Bq <C JSperc for ^ ^s, in practice their values are 
very close, since Bo/-Bperc — (10^s/^l)^/^. The predicted 
negative MR may be used for the experimental deter- 
mination of the ratio Il/Is- Strong negative MR has 
been observed in very-high- mobility samples |jl^,|l^, in 
qualitative agreement with our theory. A more detailed 
experimental check of our predictions would be desirable. 

In conclusion, we have studied the quasiclassical mag- 
netotransport of a 2DEG with smooth disorder and rare 
strong scatterers. Interplay of these two types of disor- 
der leads to a novel mechanism of strong negative MR; 
the latter is shown to saturate with increasing i? at a 
value determined solely by smooth disorder. The results 
are relevant to experiments on transport in dilute antidot 
arrays, as well as in high-mobility heterostructures with 
background impurities or interface imperfections. 

We thank M. Heiblum, J. Smet, V. Umansky, and 
D. Weiss for discussions of experimental results, which 
motivated this research. This work was supported by 
the SFB 195 and the Schwerpunktprogramm "Quanten- 
Hall-Systeme" der Deutschen Forschungsgemeinschaft. 



* Also at Petersburg Nuclear Physics Institute, 188350 

St. Petersburg, Russia, 
t Also at A.F. loflFe Physico- Technical Institute, 194021 

St. Petersburg, Russia. 
[1] C.W.J. Beenakker and H. van Houten, Solid State Phys. 

44, 1 (1991); D.K. Ferry and S.M. Goodnick, Transport 

in Nana structures (Cambridge University Press, 1997). 
[2] A.D. Mirlin et ai, Phys. Rev. Lett. 83, 2801 (1999). 
[3] M.M. Fogler et ai, Phys. Rev. B 56, 6823 (1997). 
[4] A.D. Mirhn, D.G. Polyakov, and P.Wolfie, Phys. Rev. 

Lett. 80, 2429 (1998); F. Evers et ai, Phys. Rev. B 60, 

8951 (1999). 

[5] E.M. Baskin, L.I. MagariU, M.V. Entin, Sov. Phys. JETP 
48, 365 (1978); E.M. Baskin and M.V. Entin, Physica B 
249-251, 805 (1998). 
[6] A.V. Bobylev et ai, Phys. Rev. Lett. 75, 197 (1995); J. 

Stat. Phys. 87, 1205 (1997). 
[7] G.M. Gusev et ai, J. Phys. Cond. Matter 6, 73 (1994). 
[8] K. Tsukagoshi et al, Phys. Rev. B 52, 8344 (1995). 
[9] G. Lutjering, Ph. D. Thesis (Stuttgart, 1996). 
[10] G. Nachtwei et al, Phys. Rev. B 55, 6731 (1997). 
[11] P.T. Coleridge, Phys. Rev. B 44, 3793 (1991). 
[12] T. Saku, Y. Horikoshi, and Y. Tokura, Jpn. J. Appl. 

Phys. 35, 34 (1996). 
[13] V. Umansky, R. de Picciotto, and M. Heiblum, Appl. 

Phys. Lett. 71, 683 (1997). 
[14] J.M.J, van Leeuwen and A. Weijland, Physica 36, 457 
(1967); A. Weijland and J.M.J, van Leeuwen, Physica 
38, 35 (1968). 

[15] E.H. Hauge, in Transport Phenomena, ed. by G. Kir- 
czenow and J. Marro (Springer, Berlin, 1974), p. 337. 

[16] D.G. Polyakov, F Evers, A.D. Mirlin, and P. Wolfie, in 
preparation. 

[17] J. Smet, private communication; M. Heiblum and 
V. Umansky, private communication. 



4 



